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Abstract. Bipolar molecular outflows have been observed and studied extensively in the past, 
but some recent observations of periodic variations in maser intensity pose new challenges. Even 
quasi-periodic maser flares have been observed and reported in the literature. Motivated by 
these data, we have tried to study situations in binary systems with specific attention to the 
two observed features, i.e., the bipolar flows and the variabilities in the maser intensity. We 
have studied the evolution of spherically symmetric wind from one of the bodies in the binary 
system, in the plane of the binary. Our approach includes the analytical study of rotating flows 
with numerical computation of streamlines of fluid particles using PLUTO code. We present the 
results of our findings assuming simple configurations, and discuss the implications. 
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1. Introduction 

Bipolar outflows are ubiquitious in nature and presumably thought to be associated with 
star forming regions in molecular clouds. Since the discovery of these bipolar outflows 
(Snell et al. (1980)) various attempts have been made to understand the physical nature 
of these phenomena. Such outflows are generally believed to occur around young stellar 
objects (YSOs) and are thought to have intimate relationship with the process of star 
formation and early stage of stellar evolution (see reviews Snell (1983), Lada (1985), 
Bachiller (1996)). Unanimous view about the process of star formation and its early evo- 
lution is yet to emerge which in turn, inevitably, makes it difficult to have a clear under- 
standing of the nature of these bipolar outflows. Observations of periodic/quasi-periodic 
variations in maser intensity from such regions add further complications (Goedhart, 
Gaylard & van der Walt (2004), Goedhart et al. (2005), Goedhart et al. (2009), van der 
Walt, Goedhart & Gaylard (2009), Szymczak et al. (2011), Araya et al. (2010)). 

Periodicities observed in the maser light curves are one of the most challenging, poorly 
understood, features. In the present work, we have tried to demonstrate that a binary 
system, consisting of a star and another gravitating object, could be a potential candidate 
to explain together both of the following observed features of the maser sources, the 
/\ . bipolarity and the periodicity in the intensity variations. 

2. The Model and Simulations 

Consider a binary system consisting of two bodies, S and P, which are rotating around 
their common center of mass, O, in a plane. Let XYZbe the inertial (fixed) coordinate 
frame in which the two bodies lie in the XY — plane with angular velocity vector in 
Z— direction. Let XYZbe the rotating (corotating) coordinate frame which rotates with 
an angular velocity same as that of the two bodies in the binary system and therefore 
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Figure 1. The inertial (XY) and rotating (XY) coordinate frames are shown in the plane of 
a binary. The center of mass of the two bodies labelled as S and P is the origin of both the 
coordinate frames and is denoted by O. 

both the bodies appear to be at rest in this frame. The origins of both the coordinate 
frames coincide and are taken to be at the center of mass, O, of the binary system. The 
units of various quantities are chosen such that the properties of the system depend only 
on a single parameter. Let the total mass (M) of the primaries (S and P) be the unit 
of mass; the distance between them (V) be the unit of distance; and the unit of time be 
chosen in such a way that the angular speed of the primaries, denoted by Q, be unity. 
Let £ be the mass of P, thus mass of S is (1 — £). In the rotating reference frame, with 
positive X in the direction of the body P, the coordinates of P and S will be (1 — £, 0) 
and (— £, 0) respectively. If (X, Y, 0) be the coordinate of an arbitrary point A, then from 
Figure 1, 



rsA = \rs A \ = V(X + 2 + Y 2 ; ?pa = \r PA \ = y/(X - (1 - £)) 2 +F 2 (2.1) 



(2.2) 



The gravitational potential in the corotating frame at the point A may be written as, 

* = _(!zi)_J_ 

Let us assume that one of the bodies, say 5, has a spherically symmetric wind very near 
to its upper atmosphere, whereas the other body, P, is interacting only gravitationally. 
If v(X,t) be the fluid velocity of the wind in the rotating frame then we may write the 
Euler equations in rotating frame for steady flow, with p and p as the fluid pressure and 
density respectively, as, 

Vp 
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where (X,r) = (X,Y,Z,r) and £1 (= ez, which is the unit vector along Z 
angular velocity of the corotating frame relative to the inertial frame. 
It can be shown that, 

(vV)B = 



Z) is the 



(2.4) 



where 



B 



1 2 

2 V 



[dp 
J J 



$ 



eff 



$eff 



$ 



2 1 ' 



(2.5) 



Therefore the quantity, S, is constant along a particular streamline for steady flows, 
although it could be a different constant for different streamlines. Noting the fact that 
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the particle paths and streamlines are the same for steady flows, we can see that B 
remains the same for a particular fluid element as it moves along a particular streamline. 
Thus, for a particular streamline, we may write 

* 2 , ( 7 



B = 9 V + ( ~~ 1 ) RT + ^ eff = constant = C ( 2 ' 6 ) 

where we have used adiabatic equation of state and note that the term J dp/ p appearing 
in equation 12.51 may be replaced by specific enthalpy (w) for isentropic evolution of fluid 
element. 7 is ratio of specific heats at constant pressure and constant volume, R is the gas 
constant and T is the temperature. As the terms v 2 /2 and {7/(7 — 1)} RT in equation 
12.61 cannot be negative, we infer from equation 12.61 that the motion of a fluid element, 
and hence the corresponding streamline, is restricted to the region where $eff < C. 

2.1. Numerical simulations using PLUTO code 

To study how the spherically symmetric wind from the body S flows in the presence 
of another gravitating body P in a binary system, we use PLUTO code. The details 
of the code may be found in Mignone et al. (2007)| (and references therein, and code 



at |http ://plutocode.ph .unito.it/[ ). We use the hydrodynamic module of the code and 
solve the equations in three-dimensional cartesian geometry. We adapt the code in the 
corotating frame of the binary in which the two bodies, S and P, appear to be at rest, 
by adding the necessary body-forces, namely, coriolis and centrifugal, to the equation of 
motion. To understand the evolution of the wind in the plane of the binary, we start with 
a high-pressure circularly symmetric region at the location of S and the fluid pressure 
outside this region being very small, the fluid particles experience force which is pointed 
radially outward from the location of S. We have performed simulations for various 
initial conditions (by chosing different values for high/low pressure/density regions) and 
for different values of the parameter £. Also, the flow need not be steady. Our simulations 
are done in two ways: (a) we set the initial conditions and study the flow. In this case the 
matter eventually flows out of the computational domain as there is no supply of matter 
from the location of the body *S, and the code stops; (b) Having set the initial conditions 
we may supply the matter at some arbitrary time intervals. In this case, as the matter is 
not completely depleted out of the computational domain, the code runs for longer time. 
How quickly the matter is depleted out of the domain depends on the initial condition. 
Thus, we may study the average properties of the wind flow (i.e. density, pressure etc.) 
by plotting density /pressure maps at different times. 

2.2. Results from PLUTO simulations and discussion 

We choose to present the result by demonstrating the evolution of density maps with 
time as seen in the corotating frame which is the rest frame of the binary system. This 
choice of plotting the density map seems relevant, as ultimately we will be interested 
in knowing the distribution of matter in space and its evolution in time to identify the 
regions which could potentially be maser emitting spots, particularly due to relatively 
high concentration of matter. From various panels in Figure 2, we see that the isotropy 
of the wind is broken very near to the binary system, as desired for density modulation 
seen for a fixed line-of-sight of an inertial observer. Although it is known that the orbits 
of a test-body are chaotic in the binary system, the problem known as the restricted 
three body problem (Poincare (1890)), it is remarkable to note that the third test-body 
being replaced by the wind (continuous matter) evolves in a similar fashion for various 
initial conditions and for different values of £ as seen from simulations (e. g., the average 
property of the wind, say, density, seems to evolve in a definitive way). Figure 2 is the 
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Figure 2. Snap-shots of density maps, as seen from the corotating frame (for £ = 0.5 and 
with arbitrary supply of matter at S), where the bodies S and P are at (—0.5,0) and (0.5,0), 
respectively. Time increases from upper- left (initial time) to lower-right. 

result of one of many simulations performed to study this problem. The anisotropy of 
the wind may be understood to be due to the shapes of the isocontours of the effective 
potential (4> e ff) which has five Lagrange-points near to the binary, and these isocontours 
tend to become Keplerian far away from the binary system. Hence one may expect that 
the outflowing matter, which need not escape, settles into Keplerian orbits depending on 
their initial velocities. These results prompt us to imagine that on an average, much of the 
matter spirals outward with certain pitch-angle, which depends on the intial conditions, 
and tends to settle in Keplerian orbits far away from the binary system, forming a torus 
skirting the spiral pattern. 

2.3. Results of variability (light- curve) simulation 

We simulate the situation discussed at the end of last subsection and present the results 
by showing light-curves of maser intensities from different locations as seen along a fixed 
sight-line of an inertial observer. As it should, the modulations in the maser intensity 
will depend on the inclination of the binary system with respect to the sky-plane, and 
so we show our results for an arbitrarily chosen inclination angle. Columns of negligible 
gradients in the sight-line velocity component with adequate concentration of matter are 
the most preferred sites of maser emission. 



3. Summary and Conclusions 

We have investigated the flow of the wind from one of the bodies in a binary system, 
and tried to understand the plausible mechanisms for modulations in the maser intensity. 
Anisotropics seen in the rotating frame will have the desired character of bipolar flows and 
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Figure 3. Light-curves (shown as line plots corresponding to the locations marked in the color 
panel, and over two orbital cycles) that an inertial observer would see, as monitored from our 
simulation of maser emission from flows in binary systems. 

also will appear to be periodic (with binary-period) to an inertial observer for relevant 
sight-lines. Further, in our model the potential maser spots in the sky-plane do not move 
and the minima in intensity variation cycle repeat at regular intervals of the orbital 
period. These characteristics are naturally produced in our model wherein the variabilities 
observed in the maser intensity is due to the density modulation resulting from the flow 
in the binary as seen from the simulations. 
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